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We apply a recently developed 1/(N — 1) expansion to the lull counting statistics lor the A^'-lold 
degenerate Anderson impurity model in the Kondo regime. This approach is based on the pertur- 
bation theory in the Coulomb interaction U, and is dilFerent from the conventional large N theories 
such as the usual 1/N expansion and non-crossing approximation. We have confirmed that the 
calculations carried out up to order 1/{N — 1)^ agree closely with those ol the numerical renormal- 
ization group at A' = 4, where the degeneracy is still not so large. This ensures the applicability 
ol our approach lor A > 4. We present the results ol the cumulants ol the probability distribution 
lunction lor nonequilibrium current through a quantum dot in the particle-hole symmetric case. 
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I. INTRODUCTION 



II. MODEL AND FORMULATION 



The universal Kondo behavior in the quantum dots 
has been studied intensively for nonequilibrium current 
[ll-Q and shot noise [7l-[l3j . Effects of orbital degeneracy 
on the universality have also been a subject of current 
interest ^i^. 

Essential features of the low-energy properties of the 
orbital systems may be deduced from the SU(Af) An- 
derson model, which for N = 2 corresponds to the spin 
degeneracy. The exact numerical renormalization group 
(NRG) approach [l^ has been applied to this model for 
small degeneracies Ai" < 4 [2^. For A^ > 4, the conven- 
tional large N theories, such as the 1/N expansion and 
the non-crossing approximation [2ll [23 | are applicable to 
a wide energy scale except for the low-energy Fermi-liquid 
regime. Therefore, at low energies alternative approaches 
are needed to explore the universal behavior for A^ > 4. 

Recently, we proposed a different type of large A^ the- 
ory, using a scaling that takes u ~ (N — 1)[/ as an inde- 
pendent variable instead of the bare Coulomb interaction 
U [2^, [ll] . The factor A^ — 1 represents the number of 
interacting orbitals, excluding the one prohibited by the 
Pauli principle. With this scaling the perturbation series 
in U can be reorganized as an expansion in powers of 
1 / (A' — 1). To leading order in 1/{N — 1) it describes the 
Hartree-Fock random phase approximation (HF-RPA). 
The higher-order corrections systematically describe the 
fluctuations beyond the HF-RPA. As the unperturbed 
Hamiltonian includes the tunneling matrix element be- 
tween the impurity and conduction bands, this approach 
naturally describes the Fermi-liquid state. 

In this report, we discuss dependence of the renormal- 
ized parameters on the degeneracy A^, carrying out the 
calculations up to order 1/(A^ — 1)^. Furthermore, we 
apply this approach to the full counting statistics [25l - 
[27| for the nonequilibrium current distribution through 
quantum dots in the Kondo regime. 



We consider the SU(Af) impurity Anderson model with 
a finite interaction U, connected to two leads (i^ — L, R): 
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N 

+ Y '"•^ {dln'^um + i^lmd,n) > (1) 
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Here, = e^^ + [N — l)U/2, d}^ is the creation op- 
erator for an electron with energy and orbital m 
(= 1, 2, • • • , A^) in the impurity site, and n^^^ = dl^^d^^. 
The operator c|j,„j for a conduction electron in the lead ly 
is normalized as {c^^.^, cl,^,^,} = 6^^/ i5mm"5(e - e'), and 
i/Vm = I-od'^VP'^ei^rn- The hybridization energy scale is 
given by A = F^ + r_R, with F^ irpvl and p = 1/(21?). 

We use the imaginary-frequency Green's function that 
is given, for \uj\ ^ D, by 



Giiuj) 



icj — + iAsgncj — Yj{iuj) 



(3) 



Here, 'S{iuj) is the self-energy due to Hu- The ground- 
state average of the local charge, {udm) = <5/7r, can be 
deduced from the phase shift 5 = eot^^(i?^/A) with 
E*j^ = + 5](0). The renormalized parameters are de- 
fined by 1/z = 1 — dT,{iuj) / d{iuj)\^=Q, = z E*^, and 
A = zA. Furthermore, the enhancement factor for the 
spin and charge susceptibilities, Xs = Xmm ~ Xmm> ^^'^ 
Xc = Xmm + (A^ — ^)Xmm'^ ^au be expressed in terms 



of z and the vertex function F 
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FIG. 1: The leading order diagrams in the 1/{N — 1) expan- 
sion. The wavy and sohd hnes indicate U and Go. The double 
wavy line represents Wbub(i'^) given in Eq. Q. 





FIG. 2: The order 1/{N — 1) vertex and self-energy. 



for m:^m' [21, [2£ 
1 



^r_„,„,„(0,0;0,0), (4) 



and U = ^^r^„';„'m(0'0;0, 0) for m ^ m' represents 
the residual interaction between the quasi-particles. 

We introduce a scaling for the bare and renormalized 
interactions, by multiplying a factor iV — 1 [23i. [23|. 
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With these parameters, the Wilson ratio R=zxs and 
that for the charge sector can be expressed in the form. 



R = 1 



N - 1 



sin^J, 



zxc = ^ — 9 s'm^S. (6) 



One of the merit of this scaling is that the perturbation 
expansion with respect to Hu can be classified according 
to the power of l/(iV — 1), taking g as an independent 
variable. In the present report we consider the particle- 
hole symmetric case: Cd = —{N — l)U/2 and S = 7r/2. 



III. EXPANSION UP TO ORDER 1/(N - if 

At the zero order with respect to 1/{N — 1), the limit 
of — >■ oo is taken at fixed g, and it gives the Hartree- 
Fock results. Specifically, at half-filling the unperturbed 
Green's function is given by Go{iLu) = [iuj + iAsgnoj] ^ 
as ^rf = 0. 

The leading order corrections in the 1/(A — 1) expan- 
sion arise form a series of the bubble diag rams, Z4]^^\^(^iuj^ ^ 
of the RPA type shown in Fig. [TJ 



fe=i 



(7) 



Here, Xo(*^) = ^ /"5F^o(*'^ + i^')Go{ii^j'), and the co- 
efficient Ak = {N- I f ^^^0 [-V(^ - 1)]^ arises from 



the summations over the orbital indices for a series of 
k fermion loops. Correspondingly, the order 1/(A— 1) 



contributions of F 



mm :m m 



(0, 0; 0, 0) arise from the first 




FIG. 3: The order l/{N-iy dia erams for the vertex function 



mm' :m' m 



(0, 0; 0, 0) for m / m' at half-filling 



0. 




FIG. 4: The order 1/(A - 1)^ self-energy at half-filling. 



diagram in Fig. ^ which give the leading order correction 
to the renormalized coupling, 
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This correction determines the Wilson ratio R to order 
1/(A-1) through Eq. Similarly, the order 1/(A-1) 
self-energy arises from the second diagram in Fig. [51 

Higher order fluctuations beyond the RPA appear first 
through the next-leading-order contributions. Figures [31 
and m show the order 1/(A— 1)^ diagrams for the vertex 
function and self-energy in the particle-hole symmetric 
case, respectively. These contributions and the higher or- 
der components from Eq. ([7]) determine the renormaliza- 
tion factor z and the Wilson ratio R to order 1/(A— 1)^. 
In the present work, we have calculated all these next- 
leading-order contributions. The outline of the calcula- 
tions were described in a previous paper (2^ . We have 
also carried out the NRG calculations for A = 4 to 
demonstrate the reliability of the 1/(A — 1) expansion. 

Figure [H shows the next-leading-order results for Ij and 
z, plotted as a function of g. We see the very close agree- 
ment between the NRG and the next-leading-order re- 
sults in the 1/{N— 1) expansion for A = 4, especially for 
the renormalized coupling g. The two curves for g, for 
A = 4, almost overlap each other over the whole range of 
g although the next-leading-order results (green dashed 
line) are slightly smaller than those of the NRG (solid 
circles). As A increases, g converges rapidly to the RPA 
value p — > ,9/(1 + g), which is asymptotically exact in 
the limit of A — > 00. We also see that the value that g 
can take is bounded in a very narrow region between the 
curve for A = 4 and that for the A — >■ c» limit. 

The order 1/(A — 1)^ results for the renormalization 
factor z, shown in Fig. [51 also agree well with the NRG 
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FIG. 5: (Color online) Renormalized parameters gt and z 
in the particle-hole symmetric case td = —{N — l)U/2 are 
plotted as a function of the scaled Coulomb interaction g for 
N = 2 (Bethe ansatz [lOl), and for A'^ = 4, 6, 8 (next-leading 
order in the 1/{N - 1) expansion). For TV = 4, the NRG 
results [l^.[3ll] are also shown with the solid circle (•). In the 
N ^ oo limit, g and z approach 5 ^ + g) and z 1.0. 



results for = 4 at g < 3.0, or cquivalently in the region 
where g < 0.8. This indicates that the next-leading- 
order results also reproduce the Kondo energy scale, A = 
zA, properly from the weak-coupling to intermediate- 
coupling regions where g is still not very close to 1.0, the 
value in the strong coupling limit. 



IV. FULL COUNTING STATISTICS 

The 1/{N — 1) expansion can be applied to nonequilib- 
rium transport at finite bias voltages V. To be specific, 
we choose the lead-dot couplings and chemical potentials 
in the leads to be symmetric = F^ and — ~M-R 
(= eV/2), and consider a steady state in the particle- hole 
symmetric case. In this case, an asymptotically exact re- 
sult was obtained for the retarded Green's function at low 
energies up to terms of order w^, T^, and (eV)'^ [1, [11], 



Lu + iA + i 



2{N-1)A 



(9) 



The average and fluctuations of the steady current at low 
energies can be described by the local Fermi-liquid the- 
ory, which follows from this form of the Green's function. 
The Fermi-liquid theory, or the related renormalized per- 
turbation theory (RPT) [s^, can also be applied to 
the full counting statistics for nonequilibrium transport. 



Specifically, we consider the probability distribution P{q) 
of the transferred charge q = (gi,(3'2' ' ' ' '^at) from the 
left lead to the orbital m of the Anderson impurity during 
a time interval T. The generating function for this prob- 



, Ajy), and can 



ability distribution is defined by x (•^) = J2q e.^^'''P{q) 
with the counting fields A = A2, • • • , Ajy 
be expressed in the form (25j . 

X(A) = (Tce^V^~ijjit[H^it)+nu{t)\ 



(10) 



Here, Tc is the time ordering operator along the Keldysh 
contour C, and the time evolution is defined with respect 
to an extended Hamiltonian 'H^(i) = "Hg + 'H^(i) + H^, 



H.c. 

(11) 



Here, the counting field depends on the path such that 
Am(^T) — foi' the forward (i_) and backward (t+) 

paths, respectively, and is switched on during a time in- 
terval [i : ^ r ^0]. 

We obtained an explicit expression of x{^)^ which is 
asymptotically exact at low energies u p to order {eVY for 
general and U , using the RPT (TtI. Il8l|. The cumulant 
for the full current can be derived from xC-^), choosing 
the counting fields to be A„j = A for all m, and then 



taking a derivative C„ = (- 



^lnx(A): 
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J u 

e 



Sir, 



12 



(1 + 2"+i) g^ 




(12) 



Here, J„ = Ne'^V/ {2TTh) to the linear- response current 
of the unitary limit, and Snn' is the Kronecker's delta. 
Equation (|12p shows that the nonequilibrium properties 
are also characterized by the two renormalized parame- 
ters, 5 and the Kondo energy scale A, at low energies 
in the particle- hole symmetric case. For n = 1, the cu- 
mulant eCi/T expresses the steady current. The uni- 
versal properties of the cumulants for n > 2 can be 
extracted from the Fano-factor inspired ratio (FFIR), 
which is normalized with respect to the backscattering 
current Jt = Ju — eCi/T [l^ : 



CP 



1 



(1 + 2"+i) 52 
N -1 



1 + 



T e ' 



(13) 
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Here, can also be regarded as the Poisson value of 
the cumulant, and the FFIR for 71 = 2 corresponds to 
the Fano factor [T6j . 

Figure [6] shows the ratio C„/C,f for the noise n = 2 
(top), the skewness n = 3 (middle), and sharpness n = A 
(bottom) as a function of the scaled Coulomb interac- 
tion g for several value of the orbital degeneracy N. 
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We see the close agreement between the NRG and the 
next-leading-order results for N = A again. As the next- 
leading-ordcr 1/{N — 1) results for g are numerically al- 
most exact for > 4, the results shown in Fig. IH] cap- 
ture orbital effects correctly. In the strong coupling limit 
g -> oo, the FFIR approach C„/C,f ^ (iV+2"+i)/(A+4) 
as the renormalized coupling converges to g — > 1 . As A 
increases, g converges rapidly to the value, g ~ 5/(1 + .g), 
for the large A^ limit, as mentioned in the above. Thus, 
for N > 8, the A^ dependence is determined essentially 
by the factor 1/{N — 1), appearing explicitly in Eq. (|13p . 
In the limit of A^ — > oo, the cumulant approaches the 
Poisson value, C„/C^ 1, as the fluctuations due to 
eletctron correlations are suppressed and the mean-field 
theory becomes asymptotically exact. 

In summary, carrying out the 1/{N — 1) expansion up 
to order 1/(A^ — 1)^, wc have obtained the cumulants 
for the probability distribution of the current through 
a quantum dot with the orbital degeneracy A^ > 4 at 
low energies. The results of the renormalized coupling 5 
show the excellent agreement with the exact NRG results 
already at A^ = 4 in the particle-hole symmetric case. 
This enable us to obtain almost exact numerical results 
for the Fano-factor inspired ratio C^^/C^ for A^ > 4. 
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FIG. 6: (Color online) The FFIR C^^/C^ for n = 2 (top), 
n = 3 (middle), and n=4 (bottom) are plotted as a function 
of g for iV = 2 (Bethe ansatz), and N = 4,6, 8, 10, 20 (next- 
leading-order results in the 1/{N — 1) expansion). The NRG 
results (solid circle) are also shown for AT = 4. 
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